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Abstract 

We consider the 4d compact U{\) gauge theory with extended action 

S = —(3 cos Op — ^ cos 26p 
p p 

We give a fuh characterization of the phase diagram of this model using the notion 
of flux. The relation with the usual monopole picture is discussed. In analogy with 
the XY model we consider the helicity modulus ^ for this theory, and show that 
it is an order parameter. Analyzing the finite-size effects of the helicity modulus 
we conclude that the transition is first-order. The value of this order parameter is 
related to the renormalized coupling Pr. We measure at the transition point and 
give a counterexample to its conjectured universal value [2]. 
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1 Introduction 



A phase transition in an infinite system is rigorously defined by a non-analyticity 
of the free energy as a function of the driving parameter (hke the temperature, an 
external field or some coupling of the theory). In practice, a useful description can be 
provided by identifying an observable (the so-called order parameter) whose vacuum 
expectation value is zero in one phase and non-zero in the other. This behavior is 
manifestly non-analytic and can sometimes be used to relate the phase transition 
with the spontaneous breaking of some symmetry (implicitly or explicitly defined) 
of the model, thus providing an appealing physical picture of the transition itself. 

The formulation of an order parameter is a highly non-trivial problem, and closely 
depends on the specific features of the theory considered. A broad class of theories 
(for a review see [2|), which includes the 4d compact Abelian theory we are going to 
consider, admits a characterization of the transition in terms of topological defects, 
which is a general name used to indicate some collective degrees of freedom with 
non-trivial homotopy group in terms of which one tries to reformulate the partition 
function of the system: the picture which emerges is a theory of interacting defects 
(sometimes one says that there is a gas of defects), which are very dilute in one 
phase and are condensed in the other. The nature of such a description is truly 
non-pert urbative. 

An open question is the following: is the choice of the excitation responsible 
for the phase transition unique? or is there some freedom in the identification of 
the topological structure involved? Here, we consider this question for the compact 
U{1) system on a hypertorus. 

In 3d, a convenient choice of defect is the pointlike magnetic monopole f-^; in 
Ad we naturally have, instead of pointlike objects, monopole currents representing 
monopoles evolving in space-time. The partition function of the theory can be 
expressed in terms of these currents jH] and the following picture emerges: in addition 
to a weak-coupling Coulomb phase in which only small and dilute loops are present, 
one finds a strong-coupling confined phase in which there are many monopole loops 
deeply entangled. 

Suppose now that one monitors the value of the electromagnetic (e.m.) fiux 
through some plane; we expect to observe two completely different situations in the 
two regimes: the flux should be constant if no monopoles are present, while it should 
dramatically vary in the other phase. This is the physical picture we are going to 
consider, namely a picture of the phase transition in terms of fluxes and not in terms 
of monopoles: they are equivalent from the logical point of view; there can then be 
arguments in favor of one or the other picture, but none 'precedes' the other in a 
strict logical sense. 

As we will extensively see, this alternative point of view provides a way to define 
a natural order parameter for the system, the so-called helicity modulus, defined by 
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the response of the system to a variation of the flux. An accurate determination of 
this observable will be used as the basis for a finite-size effect analysis, which will 
provide good evidence about the first-order nature of the transition. 

The paper is organized in the following way. In sections |21 and 0] we review the 
model under study and the notion of flux on a lattice. In section 0] we show how the 
behavior of fluxes qualitatively changes across the phase boundaries of the system, 
and thus provide a first characterization of the model in terms of fluxes. We also 
interpret this behavior in terms of the monopole picture, showing the equivalence 
of the two approaches (we will not review the monopole picture |3] j3] , but only use 
some results to interpret our data). In section El we develop an analogy between 
our Abelian context and the non-Abelian one, proposing the analog of 't Hooft's 
twisted boundary conditions jB]. In section IHl we introduce an order parameter for 
our theory, the helicity modulus [7j. In section [7| we show our numerical results for 
the order parameter in the whole phase diagram; in the following three sub-sections 
we discuss {{) the issue of the order of the phase transition, (ii) the relation between 
the helicity modulus and the renormalized coupling (in the Coulomb phase), and 
[iii) the conjecture (first proposed by Cardy [2J) about the universahty of its value 
at the transition. In section |H1 we show how the presence of a monopole in the 
lattice can be enforced by adding a flux and using C-periodic boundary conditions; 
the free energy of a monopole is studied and a comparison with other constructions 
is presented. Conclusions follow. Preliminary results about sections 2-5 have been 
presented in jH]. In TableElin the Appendix we present a summary of the simulations 
performed for the paper. 



Consider a 4d lattice with periodic boundary conditions (p.b.c). Call a the lattice 
spacing. We associate to links a f/(l) field U^{r) = exp(i6'^(r)) = exp{iaeAf^{r)), 
where e is the gauge coupling, r = ax^e^ {x^ = 1, . . . , L^) and is the unit vector 
in direction /z = 1, . . . , 4. We work with the following action: 



where 9p = 6^^{r) = 6^{r) + 6y{r + e^) — 6^{r + e^) — 0^{r) is the plaquette angle, j3 
is the Wilson coupling, and we refer to 7 as the extended coupling. This model was 
introduced by Bhanot [9^ and has attracted interest also recently |10j. 

In order to point out the symmetries of the model consider the following trans- 
formation of link variables 



2 The model 




(1) 
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Ui{x,y,z,t) 
U2{x,y,z,t) 

U3{x,y,z,t) 



Ui{x,y,z,t) 

U ■ U2{x, y, z, t) if X is odd 

U2{x,y,z,t) if X is even 

U ■ U^{x, y, z,t) if X + y is odd 

U^{x,y,z,t) if X + ?/ is even 

tf • U^i^x, y,z,t) if X + y + z is odd 

(/^{x, y,z,t) if X + y + z is even 



Ui{x,y,z,t) 

where f/ is a U{1) matrix at will. This transformation is such that 



(2) 



• if [7 = — 1 = e*'^ every configuration with couplings 7) is mapped into 
one with couplings (— /5, 7). But since this is just a change of variables, we 
conclude that the mapping (3 — > —(3 is a symmetry of the theory. 

• ifU = we map (0,7) — > (0, —7), which is therefore another symmetry. 

Furthermore, the change of variables d^{x) 29^{x) maps the couplings (/3,0) onto 
(0,/3), showing the equivalence of the model on the two coordinate axes. 
Finally, the limit 7 — >■ +00 corresponds to recovering the Z2 gauge theory {6p = 
{0, vr}; a gauge-transformation then fixes the links to Z2). 

In Fig. ^we draw a sketch of the phase diagram. 
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Figure 1: Phase diagram of the extended U{1) model. 
Three unknown aspects of the phase diagram are indicated with question marks: 

3 



1. the nature of the phase indicated with //; 



2. the nature of the phase IV; 

3. the fate of the two phase boundaries in the bottom right corner of the phase 
diagram (the same holds for the bottom left corner). 

Another long standing question, still not completely resolved, is the order of 
the phase transition, especially along the / — /// line. A recent paper ^T] shows, 
through a high statistics analysis of the plaquette action distribution, that the tran- 
sition is first-order along the Wilson axis. We will confirm this result using a totally 
different observable. It is also numerically known that the transition becomes weaker 
and weaker as 7 becomes negative, which motivated the hypothesis that a tricritical 
point could exist ^2] along the phase boundary between regions I and /// (see 
Fig. Q), beyond which the transition could be second-order. We will also address 
this issue with our new observable. 

We now introduce the notion of flux, which will provide the tool we need to 
explore the phase diagram. 

3 The flux in Abelian gauge theories 

To introduce the notion of flux we need in Ad, it is more natural to start from the 
2d case. Consider a compact U{1) gauge theory defined on a 2d lattice, with sizes 
L^,Ly, and p.b.c. Let us consider the definition 



where [6p\^.„^t, is the plaquette angle reduced to the interval (— 7r,7r) (its so-called 
physical part), and the sum is over all plaquettes. This quantity, due to p.b.c, is 
2T[k valued {k & Z), because each link is summed twice with opposite signs, giving, 
overall, modulo 2tt. 

If we consider the limit a ^ 0, we obtain the standard result 




(3) 



p 



(4) 



where F^^, is the field strength; so the definition in Eq. ^ becomes 




(5) 
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which motivates the interpretation of $ as an e.m. flux. Furthermore, in this hmit 
we notice that fluxes corresponding to different values of k are topologically discon- 
nected, in the sense that no tunneling is possible between two flux values through 
a smooth variation of gauge field. In fact, to change the value of the flux, one pla- 
quette angle must exceed the value of vr; the statistical suppression of this change is 
Q-i^^s e~^^, which goes to zero in the continuum limit (/3 — oo). ^ 

Consider now the same theory, but with (1 = 4. Call L^, L^,, L„ and Lp the four 
sizes of the lattice. Our definition of flux through any v) orientation is: 

^ " iiv planes 

A double sum is present: 

• the internal Y2p is the sum over the plaquettes in a single plane. So plane 
by plane we use the same definition as in the 2d case. 

• the external average pianes parallel planes of a given orien- 
tation is non-trivial because the flux through different planes can change (in 
3d, for example, there are monopoles between planes). The allowed values for 

are thus multiples of 2'7i/LpL„. 

Configurations whose flux is 27ck play a special role, for reasons which will be 
clear in the next section, and we say that they define flux sectors. 

Let us compute the action of a flux sector in the classical (minimal action) limit. 
In this case the flux through a single plane is equally distributed over all plaquettes 
and it does not change across parallel planes: 



keZ (7) 



{yp)f.u - -j—j- - YY 

We obtain (here for the Wilson action) 

S, = 2Pn'k'^ (8) 

showing that higher values of the flux are strongly suppressed. We stress that in 
the classical limit, even if $/27r is not integer valued because of the presence, e.g., 
of an external field, the action is anyway quadratic in the flux: 

S = (9) 

^In the 2d case the lack of ergodicity through flux values of simulations at large /3 can be cured 
by a cluster algorithm 
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Figure 2: Flux effective potential in the Coulomb phase at j3 = 1.2 and 7 = (for 
lattice size L = 4,6,8). The dashed parabolic line is Eq. © with (3 = 0.85; for a 
discussion of this value see section 17.21 



4 Characterization of the phase transition 

In this section we consider the behavior of flux sectors across a phase boundary. For 
this purpose we measure the flux distribution p{(j)) through one (/z, z/) orientation 
during a Monte Carlo simulation. Since this quantity extends over several orders of 
magnitude, we will show the flux effective potential w{(j)), defined by 

p(0) = e-'"^'^) (10) 
The results of our simulations are the following: 

1. In the Coulomb phase we observe minima of the flux effective potential at flux 
values 27rfc [k G Z) (Fig. As the thermodynamic limit is approached we 
observe barriers of diverging height between those minima, while the flux ef- 
fective potential at 27ik does not change. Extrapolating to the thermodynamic 
limit, we obtain that no tunneling is allowed between sectors: this property 
finally motivates the use of the word 'sector' we introduced at the beginning. 
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Figure 3: Flux effective potential in the confined phase aX (3 = 0.8 and 7 = 
(L = 4,6,8). A parabolic interpolating curve is superimposed (p(0) = e"*^ /^'^ , 
a = 2.8). 

2. In the confined phase (Fig. ^ sectors are not defined, in the sense that no 
local minima of the flux potential are observed. The distribution of the flux 
values is continuous (in the limit L oo), gaussian and independent of the 
lattice size. 

From the point of view of the flux picture of the model, Figures |21 and El characterize 
completely the phase transition. 

We now want to interpret all these results in terms of the monopole picture. 
We first discuss the interpretation of a flux tunneling event (in the Coulomb phase, 
when L < cxo, we observed that tunnelings across flux sectors are not forbidden, only 
suppressed). Let us start from the observation that monopoles are created in pairs 
and that a plane pierced by a Dirac string 0| gets an extra flux of 2tt. Consider a 
monopole configuration of the kind shown in Fig. Elwhere a monopole-antimonopole 
pair at a certain distance d wraps around the periodic boundaries; in Fig. |^ we show 
that a stack of d planes is pierced by a Dirac string. Let d — > (the lattice size): 
due to p.b.c. the two particles annihilate, and the configuration realized is the so- 
called Dirac sheet. This is the configuration responsible for a tunneling event; any 
other configuration produces a variation of fiux only in a fraction of the planes, and 
thus is responsible for the non-integer {mod 27r) values of the fiux. 

In the Coulomb phase only a few, small loops are created, so it is very unlikely 
to have monopole lines wrapping around the boundaries. If the volume is finite, 
there is a finite probability p for this event to happen, which vanishes like as 



7 



-» » 



Figure 4: Monopole configuration corresponding to a flux tunneling event. The 
dashed hne represents the Dirac string, the sohd hues the monopole world-lines, 
wrapping around the boundaries. 

L — > oo, as one can observe in Fig. |21 The reason for this exponential decay is that 
a tunneling event is possible only in presence of two monopole loops of length L 
(wrapping loops), and the free energy cost of a loop is proportional to its length. 

In the conflned phase, instead, a large number of monopoles and percolating loops 
of arbitrarily large size are present, so that the flux values through different planes 
are totally de- correlated] the flux through one plane can be considered as a random 
variable with a given distribution. The same deflnition applied to another plane 
gives another independent random variable, with the same distribution. Therefore 
the total flux Eq. (jH)) is just the normalized sum of independent random variables, 
whose distribution is gaussian as L — > oo. It is remarkable that also the smallest 
lattice we used (L = 4, that is, with 16 parallel planes in the /i, v orientation) has 
the same distribution as the largest one, indicating that the central limit 'converges' 
quickly to its asymptotic distribution. 

Moreover, the variance a of the gaussian distribution is related to the density 
5 of monopole currents. Call / the fraction of monopoles which contribute to flux 
disorder. Then, between two successive planes, L'^f6 monopoles contribute to flux 
change. If the global flux (averaged over planes) has variance cr^, the flux through 
one plane has variance L^cr^ and the difference of flux between two planes has vari- 
ance 2L^cr^. In this way we estimate 



where we divide a by 27r to relate the units of flux (g 2txN) to the number of 
monopoles (g N) through which we deflne 5. Knowing a (see caption of Fig.|2I) and 
measuring 5 = 0.488, we estimate / ~ 0.81: the physical picture is thus that 80% of 
the monopoles contribute to flux disorder at P = 0.8. We checked that this result 
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does not significantly change if we consider stronger couplings. 

We now briefly describe the algorithm we used. The measurement of p(0) 
(Eq. (fTUj) ) is almost impossible using an ordinary local updating algorithm, be- 
cause of the strong suppression of higher flux values we mentioned at the end of 
section |21 To overcome this problem, we implemented a multicanonical algorithm 
|14j : iteratively, starting from an initial distribution of flux values over a region 
Rq = {4> \ ~4>o < (f> < ^o}; we biased the sampling probability of the flux such 
that we could extend the region over which is known while maintaining a flat 
flux distribution in the middle. This bias is removed during the statistical averaging. 

Finally, in order to increase the flux tunneling probability, we implemented a 
global update of the lattice in the following way Consider the LpL„ planes in the 
orientation (/i, u), and for each plane select L^Lj, links, as shown in Fig. El Following 
the path we draw, we add to the k-th link the quantity j^j^k, k = 1, . . . , L^L^,, so 

that to each plaquette in the plane an extra flux j^yt imposed. Overall, we add 
a flux 27r through the chosen orientation. To satisfy detailed balance the orienta- 
tion of the extra flux is chosen randomly, and a Metropolis step is implemented. 
Although the acceptance probability of such a global step is normally very small, 
in our multicanonical approach it was greatly enhanced, thus providing a valuable 
tool to decorrelate the flux faster. 
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Figure 5: Picture of a plane through which an extra 27r flux is imposed. The 
set of L^Ljy links that are modified is indicated with thick lines. The quantity 
j^j^k, k = 1, . . . , L^Ly is added to the k-th. link crossed by the dashed line. 



5 Twisted boundary conditions: U(l) vs. SU(N) 

We extend to the Abelian context the idea of 't Hooft [HI about non-Abelian gauge 
theories: we probe the response of the system to a variation of the flux. In the 
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context of SU{N) pure gauge theories, we know that it is possible to modify the 
flux (by twisting the boundary conditions) according to the discrete group Z^, center 
of SU{N). This modification arises from the requirement of changing the boundary 
conditions of a gauge system defined on a hypertorus in the most general way 
compatible with gauge invariance. We can implement a similar change also in the 
Abelian context, but now the difference is that we can vary the flux continuously. 
For this reason we will intentionally use the word 'twist' as a synonym of 'fiux'.^ 

In order to implement this idea consider a stack of plaquettes, one in each plane 
of a given orientation, as indicated (for the 3d case) in Fig. IHl and change their value 
as: 

Op^ep + cf), (j)eR (12) 




Figure 6: 3d representation of the 'twists' in U{1) gauge theory. Each (yU, i^) plane 
now carries flux 0. 

The new partition function is 

Z(„ ^ / ««e=„..-"..-.^.»<-.-..^=^.--.^.»-.. (13, 

where stack indicates the set of plaquettes through which the extra flux is imposed. 
For instance 

stack = {6^i,{x, y,z,t) \ fi = l,u = 2;x = l,y = 1} (14) 

'stack^ is its complement, that is all the other unchanged plaquettes. Observe that 
Z((/)) is clearly 2n periodic, so actually the flux we add to the system is defined only 
mod 2n. It is not manifestly translationally invariant, but it is easy to show that we 
can freely move the position of the stack through a redefinition of the link variables. 
In Fig. [7|we show a two-dimensional cartoon of this change of variables: u is a link, 
and (f) inside the plaquette indicates the extra flux. A change u —>■ e~'^^u moves 
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Figure 7: Graphical representation of the change of variables that moves the twisted 
stack, u indicates a link; is the twist imposed on the plaquette. 

the stack to a neighboring position. With a similar transformation (namely, moving 
only part of the flux from the initial plaquette) one can spread the flux arbitrarily 
over several plaquettes, in particular, homogeneously through the whole plane. 

Observe that we can implement the presence of an extra flux also in terms of 
modified ("twisted") boundary conditions 

U2{x + Li, y = 1, t) = t/2(x, y = l,z, t) e"^ (15) 
This remark will play an important role in the following. 

We measure the free energy of the flux, defined by 

i^(0) = -logf|| (16) 

and the results are as follows: 

1. In the confined phase (Fig. |H1 lowest curve), within our statistical error, the 
free energy is independent of 0. The reason for this phenomenon is clear 
from Eq. (jl5|) . Since the flux variation can be interpreted as a modification 
of the boundary conditions, we expect that in the confined phase the system 
is insensitive to them, because the correlation length ^ is finite (finite size 
corrections to F{(j)) are (9(e~^/^)). 

2. In the Coulomb (massless) phase the interaction is long-range and the system 
becomes sensitive to the boundary conditions (Fig. |H1 upper curves). 

'^For an earlier, related study of compact C/(l) in an external field, see [T?)j . 
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It is possible to describe this behavior in terms of classical arguments: we 
write the partition function of the system as a sum over flux sectors of the 
classical partition functions corresponding to each sector, as derived in Eq. (jHJ 

Z{<P) = J2 Zcd<P + 27rfc) = J2 e-^(<^+^-'=)^ (17) 

k k 

We fit the flux free energy according to Eq. (fT7|) . letting /3eff be a free param- 
eter. The good description of the data provided by this ansatz (see Fig. |H1) 
motivates the interpretation of (3es as an effective coupling, which takes into 
account the quantum fluctuations (that is, the contribution of the monopole 
excitations) of the system. In section El we will discuss this point in more 
detail. 




Figure 8: Free energy F{(f)) in the Coulomb phase {(3 = 1.5, 1.1, upper and middle 
curve) and in the confined phase {(3 = 0.8, lower curve) (4"^ lattice). The fit by the 
ansatz Eq. (|T7jl is superimposed to the data. The effective couplings in the Coulomb 
phase are PeS. = 1-2, 0.71 respectively. 



3. In the unknown phases II and IV we observe (within the statistical error) the 
appearance of an extra vr periodicity (Fig. E)). We now propose an argument 
to interpret this finding. 

Consider the following Polyakov loop correlator 

{P,ix)Pi{x + L,e,))^ (18) 
12 




Figure 9: Typical free energy profile in phases II and IV {(3 = 1.0, 7 = —1.5; 4^ 
lattice). The fit by the ansatz Eq. (fTTj) (with the change 27r — tt) is superimposed 
(giving a value of /?efr ~ 3.98). 



Here Pq is the Polyakov loop of charge q in direction z> 

Pg(x) = e'«^'=i^^("+''^^) (19) 

We take the average w.r.t. Z = J^^ d(j)Z{(f)), in order to implement the influence of 
all possible external fluxes {mod 2tt) on the system; the physics anyway does not 
change. Even though the correlator is measured between point x and its periodic 
image x + L^e^, it is non-trivial due to the 'twisted' b.c. (Eq. (jl5j) ). Gauss law gives 

P,(x + L^e^) = e-*^<^P,(x) (20) 

that is the two Polyakov loops differ by the flux we added. Inserting this relation in 
Eq. (HHD we get 

{P,{x)Pl{x + Le,)), = ^VJ 7^ (21) 

Jo d(j)Z{(f)) 

If Z{(f)) has periodicity vr, it follows that when g = 1, or more generally odd, the 
numerator vanishes and the correlator is zero, while if q is even, it can be different 
from zero. Noting that 

{P,{x)P^(x + L,e,)) ~ e^^''^^'^)^'' (22) 
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where Vg is the static potential between the static charges ±g represented by the 
Polyakov loops, a vanishing Polyakov loop correlator implies an infinite (confining) 
potential. 

We therefore claim that the two unknown phases II and IV are just Coulomb 
phases for even electric charges, while odd charges are confined and paired in even 
combinations. 

Maybe a simpler argument can start from the observation that this result is 
trivially true in the Coulomb phase along the axis (3 = 0, where the vr periodicity is 
a simple consequence of the symmetry 0) (0,/5) discussed in Sec. El It then 
easily extends to the whole regions II and IV of Fig^ since a phase is characterized 
by a given infra-red (IR) behavior. Indeed our numerical results (obtained at /3 = 1.0 
and 7 = —1.5) (Fig. EI) show no detectable deviation from tt periodicity. 

The extra periodicity is also a useful tool to discuss the fate of the two phase 
boundaries at the lower corners of Fig. [TJ It is possible that these phase bound- 
aries meet at some point or that they just become asymptotically close to each 
other. What we can exclude is that they meet and terminate somewhere, letting the 
two different Coulomb phases communicate. Suppose that they communicate some- 
where, then we expect the twisted free energy F{(f)) (Fig. IHlandlHI) to be the same in 
both: but this leads to a contradiction. In fact F{(j)) would have least periodicity vr 
and 2tt at the same time. This is only possible if the distribution is fiat, which only 
happens in the confined, not the Coulomb, phase. 

We conclude this section with some comments about the numerical strategy used 
to compute the free energy Eq. ()16p. Here we are faced with an overlap problem, 
that is, in the ratio -fj^; the numerator and the denominator sample non-overlapping 
regions of the configuration space. We can cure this situation through the following 
strategy ^7]. Divide the interval [0, 27r) into parts, and consider the identity 



where is defined as the partition function of a system with an extra fiux = ^k. 
We observe that for each ratio the overlap between numerator and denominator is 
much better now. A further improvement of the overlap (and, therefore, further 
variance-reduction) is obtained if we express each ratio as 



Z{27i) 



Zn Zn-i 



(23) 



Z(0) 



Zn-i Zn-2 



Zo 



Zk _ Zk/Zi,_i 



(24) 




where the index k — ^ refers to an external fiux = (A; — |) The numerator and 
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the denominator of each ratio are ratios themselves, where, for example. 



^ = (e ^ 1 (25) 

Zju 1 

K 5 

S'fc is the action of the system when an external flux /c ^ is added. This ratio 
can be directly computed. The algorithm consists of measuring all these ratios 
separately and recombining them according to Eqs. (j^Hj) ()24|) . It is a trivial matter 
to extract statistical errors since the N ratios are estimated by independent Monte 
Carlo simulations. 



6 The helicity modulus 

We now go back to Fig.|Hl i.e., to the free energy of a system with twist 0. Consider 
the curvature of the free energy profile at some arbitrary value of the flux = 0. 
In the confined phase it is identically zero (in the thermodynamic limit), because of 
the independence of -F(0) from 0. In the Coulomb phase this quantity is different 
from zero. Therefore there is certainly a point of non-analyticity for this quantity 
at the transition: it is an order 'parameter. 
We are then motivated to introduce the following definition 

MA7) = ^^U=^=o (26) 

known as the helicity modulus in the literature PP |2| (we stress that the choice 
= is arbitrary). This name was conceived in the context of the XY (rotor) 
model, where the boundary conditions are changed by adding an angle to the 
rotor phase 0{x,y) 

e{x + L^,y) = e{x,y) + (j) (27) 

so that the quantity (j26| assumes the meaning of the response function ('modulus') 
to a rotation ('helicity') of the boundaries. Eq. ipTjl is the analog of Eq. (|TH|) in our 
case. 

If we substitute our partition function Z{(f)) Eq. (fT^ into Eq. we obtain the 
following expression 



hiP,l) = ( J](/?cos[^p] +47cos[2^p])) - (( J] (/3sin[^p] + 27 sin[2^p]))2) (28) 

stack stack 

where 'stacF refers to the stack of plaquettes we changed by adding the flux (see 
Sectional). As we noted earlier, we are allowed to spread the flux homogeneously 
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through the planes via a coordinate transformation. If we add to each plaquette of 
a given orientation the flux cf)/ L^^Ly and recompute Eq. ()26|) . we obtain equivalently 

M/?,7) = T^-^{(J](/3cos[^p]+47cos[2^p])) 

- ((5^(/?sin[^p]+27sin[2ep]))2)} (29) 
p 

where now the sum extends over all plaquettes in the (/i, v) orientation. This ex- 
pression turned out to give a much smaller variance than the other one. 

Suppose that 7 = 0. Then it is interesting to observe that in the confined phase, 
where this order parameter is zero, the average action (first term) is related to a 
measure of its fluctuations (second term). 

7 Helicity modulus: numerical results 

We now present the results of a numerical study of the helicity modulus. First, 
in Fig. we show the behavior of this quantity on the Wilson axis. Remaining 
for the moment on a qualitative footing, we see that it indeed behaves like an 
order parameter. On the extended line {[3 = 0) we would obtain the same plot 
up to a factor 4, as is clear looking at Eq. ()28|) and remembering the symmetry 
(/5,0) ^ (0,/5) of the partition function Z(/3,7). 




Figure 10: Helicity modulus on the Wilson axis near the phase transition. 
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The statistics used for this measurement is of about 10^ measurements per lat- 
tice size, with one heatbath and four overrelaxation sweeps between two subsequent 
measurements. The curves in Fig. are obtained through a reweighting of the data 
a la Ferrenberg-Swendsen [TH] . 

Let us now consider the behavior of the hehcity modulus in the whole phase 
diagram. Again refer to Fig. Q fix the value of f3 = 2.50, for example, and consider 
the (negative) values of 7 in a region across the phase boundaries of the diagram. 
We show the result in Fig. ^2 Also in this case the expectations are qualitatively 
satisfied: the fiat region signals the presence of a confining phase, and the helicity 
modulus is roughly 4 times higher in the even-charge Coulomb phase than in the 
other one. 



4 
3 
2 
1 




L=6 
L=8 
L=10 



IV 



-0.8 



-1.2 



-1.6 
Y 



-2.4 



Figure 11: Helicity modulus versus 7 (/? = 2.5 fixed; L = 6, 8, 10). Label '/' indicates 
the usual Coulomb phase, '///' the confined phase, 'IV' the even-charge Coulomb 
phase. 

From the numerical point of view, even if the statistics used is roughly the same 
as on the Wilson axis (at least across the phase transitions), these data are less 
accurate than the other ones: at such a small coupling the Monte Carlo dynamics 
is slower (much larger autocorrelation time). 

In the next sub-sections we perform a quantitative analysis of these data, in 
order to clarify the issues of the order of the phase transition, of the relation between 
the helicity modulus and the renormalized coupling, and finally of the conjectured 
universality of the helicity modulus at /3c 
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7.1 Order of the phase transition 

The difficulty of determining the order of a phase transition using lattice simulations 
is well known, and, in general terms, is due to the fact that in finite systems no 
singularities of the free energy can be present. In our specific case we are interested 
in distinguishing a ffist- from a second-order transition, but we know that even in 
the ffist-order scenario the correlation length at the transition point can be quite 
large, so that for lattices of size L < the transition mimics a second-order one. 
This can explain the difficulty that was encountered in solving this problem so far. 

As shown in Fig. ^1 the finite-size (FS) effects act to round off the singularity 
into smooth curves. It turns out that their functional form depends on the order of 
the phase transition, so that an ansatz for the FS effects, when available, provides a 
direct way to prove or disprove a conjecture about the nature of the phase transition. 
Since we have accurate data for the helicity modulus for lattices of size up to L = 14, 
we try to make such an analysis both for the ffist- and the second-order scenarios. 

From a theoretical point of view, little is known about the order of the transition. 
If the correlation length would diverge at the transition point (as happens in the 
second-order scenario), it would be possible to extract a continuum limit, leading 
to a continuum confining Abelian theory, which seems unlikely. From a numerical 
point of view, a quite recent study [n] of the id Abelian Wilson action strongly 
supports the ffist-order scenario, through a high statistics analysis of the action 
density distribution (considering lattices of size up to L = 18). 

Let us ffist consider the ffist-order scenario, and make the hypothesis that the 
helicity modulus has a jump across the transition point. We now consider the 
simplest ansatz for the partition function in the vicinity of the transition point 

^^e-^^-(^)+Xe-^^+(^) (30) 

where /+ denote the free energy densities of the two coexisting phases (such that 
f-iPc) = f+{Pc)), and X is the so-called phase asymmetry parameter, that is the 
relative weight of the two phases at the coexisting point. Eq. is the Borgs- 
Kotecky ansatz ^H] [20] , but we have to observe that the work of these two authors 
is based on a Peierls-like contour picture of the system at the transition point, with 
the hypothesis that the contours (that is, the boundary of domains characterized by 
a single phase) do not interact with each other. The long-range nature of the pho- 
ton interaction in the Coulomb phase makes it difficult to justify such a hypothesis. 
A posteriori, we will indeed observe deviations from the predictions based on this 
ansatz: for example, finite-size corrections to the ffiting parameters different from 
the standard ones proportional to ^ will be necessary. 
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The ansatz of Eq. (j3U|) leads to a formula for the interpolating FS curves which 
we write in the following way 



^(/^'^) = TTx^Ti=WTO (31) 

where is the value of the helicity modulus in the pure Coulomb phase (see Fig.lT^. 
and A/' = — f'_ is the latent heat, obtained from the expansion 

f±iP) = fc + f±-if3- Pc) + OiA/3'). (32) 

Since for simplicity /i+ is taken as constant in the Coulomb phase, the purpose 
of the ansatz is to describe the helicity modulus just below the transition, and we 
perform a fit for values of (3 < (3c (actually the range was limited to an interval 
[/?i,/?c], with Pi close to Pc, to take into account the approximation introduced by 
linearizing Eq. (j221))- For Pc we use the very precise value in [TT], and also we 
determine it self-consistently from the fit of our data (iteratively, we fit the data in 
an extended P range, determine Pc and make another fit in a range up to the new 
value; we repeat this procedure until Pc stabilizes). An analysis of the possible FS 
corrections indicates that a minimal choice is the following 



h+ ^ /i+-(l + -^) (33) 
A/' - A/'.(l + ^) (34) 

For the latent heat we consider 1/V corrections, standard for a first-order transition. 
To justify the 1/L corrections used for h^, first let us look at the behavior of the 
helicity modulus at the transition point fFig. [T^. 
We fit these data with the ansatz 

KiL) = hcioo) ■ (1 + ^) (35) 

where c is a constant, and obtain the result u = 1.04(12) , which is a first numerical 
indication for the choice of the FS correction Eq. (j33|) . At the end of this paragraph 
we are going to comment more extensively about this choice. Then, fixing z/ = 1, we 
get the thermodynamic value hc{oo) = = 0.371(1) (c = 1.7), independently 

from the Borgs-Kotecky ansatz. This value is only slightly increasing if we exclude 
from the fit the points corresponding to L = 4, 6. 

We then turn to the ansatz Eq. (pTTjl and fit the data corresponding to lattices 
L = 8, 10, 12, 14, obtaining the following results: 
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Figure 12: Helicity modulus at the transition point versus lattice size. 
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Comparing the quantities /?c, logX and A/'/6 (the latent heat per plaquette) with 
earlier measurements by Arnold et al. the agreement is very good for all 

the values. 

As a crosscheck, we can relate the value of hc{oo) = 0.371(1) measured before, 
and /i+ = 0.381(5) in Eq. (jHEl)- Their physical meaning is explained in Fig. 1131 /i+ 
is defined as 



h+ = lim^ h{p) (42) 

and its value depends on the Coulomb phase only, while he is smaller, and takes into 
account the coexisting confined phase. As one can easily check, the relation 

M/5c) = (43) 
(Eq. (j3T| computed at /3 = Pc) is consistent with our numerical values. 
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Figure 13: Graphical explanation of the physical meaning of the parameters /i+ 
and he = h{Pc, V = oo), as defined in Eq. (jHH). 

We also measured the value of h+ directly, extracting from the canonical ensem- 
ble at Pc the subset of configurations whose action density is lower than a threshold 
value (0.359) chosen to best separate the two phases; in this way one can probe the 
Coulomb phase contribution only. We find = 0.4325(15), a value incompatible 
with the value fitted using the BK ansatz. The discussion of this problem is post- 
poned to the end of this section. 

All the errors are computed performing a jackknife analysis. The values of ai, 
are mentioned, to give an idea of the magnitude of the correction involved. Since we 
want to show the relative correction that they induce, ^2 (which has dimension L^) 
is expressed as a fourth power. To evaluate the quality of the fit we could not use 
just the reduced because the data we fit are re- weighted, so that the number of 
points in the fit and their mutual correlation play a role which is not properly taken 
into account by the value of x^- So we decided to show directly a plot (FigHH left) 
of the fitted data. The most important signature of the first order nature of the 
transition is the exponential (and not power law) behavior of the helicity modulus 
as a function of the distance from the transition point, which is visible in the linear 
behavior (in log-scale) away from /3c (see Fig. IT^ . 

Now, for completeness, we want to compare the quality of this ansatz with that 
of the second-order one. If the transition would be second-order, all quantities (in 
the vicinity of the critical point) should satisfy the scaling hypothesis, that is they 
should be a function oi Lj^{^'vs, the correlation length) only. Since 

e~(/?-/?c)-'^ (44) 
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Figure 14: Graphical comparison of the quahty of the first and second order ansatz 
for the hehcity modulus data. The solid lines (left) show the first-order fit Eq.(|^. 



a rescaling of (3 according to 

f3^L'/'^{f3-f3,) (45) 

should produce a superposition ('collapse') of the data. The result of this analysis 
is shown in Fig. El on the right for the 'test' value of u = 0.3. The data collapse is 
quite poor: we consider this as an indication that the second-order scenario can be 
excluded. 

We conclude this section with a discussion of the inconsistency we found between 
the determinations of /i+ using the BK ansatz {h^ = 0.381(5)) and the direct mea- 
surement (/i+ = 0.4325(15)). On one hand, both methods are affected by systematic 
errors. Concerning the direct measurement, notice that a single metastable phase is 
not really a pure phase, but contains small 'bubbles' of the other one; furthermore 
the method is sensitive to the threshold value, which is chosen somewhat arbitrarily 
between the two peaks of the action density distribution. However, these two kinds 
of error diminish as the thermodynamic limit is approached; thus we were led to 
complement the other simulations with one on a large lattice (32^) in order to con- 
sider these errors as small. Let us now scrutinize the BK ansatz. First of all we can 
consider the possibility that the lattices we used are too small, so that significant 
higher order 1/L^,k > 1 corrections should be added to the leading terms; but if we 
exclude from the fit the smallest lattices (6^,8^, 10^), no clear increase in the fitted 
value of /i+ is observed. Similarly, no clear change in the value of the fitting param- 
eters occurs if we include a FS correction for X, the phase asymmetry parameter, 
not considered in Eqs. (j^Sl) (EH) • It appears that the BK ansatz does not describe so 
well the data very near the transition. The reason should be found, as we already 
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mentioned, in the long-range nature of the photon interaction, which is not taken 
into account by this ansatz based on a strong first-order picture. Note that correc- 
tions of order up to l/V^ were found necessary in to describe the behavior of 
the pseudo critical Pc{L) on lattices of similar sizes to ours. From this discussion we 
also obtain some plausible motivation for the unusual FS correction {1/L instead 
of 1/V at the leading order (Eq.(j2ni))) we observe for the hehcity modulus at the 
transition; we can impute such large correction to the interacting interfaces between 
the coexisting phases. 

Therefore, we blame our inconsistency between the fitted and the measured 
value of /i+ on the limitations of the BK ansatz. The direct measurement of /i+ is 
independent of any specific model of the system at the transition, and the systematic 
errors that affect its value are certainly decreasing as the thermodynamic limit is 
approached. Thus, whenever possible, we will use this determination. 

7.2 The helicity modulus and the renormaUzed couphng 

We now want to discuss a remarkable property of our order parameter |i2j: the 
helicity modulus in the Coulomb phase is related to the value of the renormalized 
coupling. Let us go back to the classical Eq. Q 



where, as already defined, F{(j)) is the flux free energy. Following [2] we can define a 
so-called 'blocking transformation'. As the number of RG steps increases and short 
distances are integrated out, only one relevant, Coulombic coupling survives, while 
all other couplings go to zero. In this way, starting from Eq. ()46|). we end up with 
the result 



where the bare coupling has been substituted by the renormalized one, but no other 
changes occurred. The second derivative of this quantity thus gives the renormalized 
coupling, as anticipated. 

The next step is to observe that the functional form of Eq. ()47|) is not totally 
consistent: this equation does not present any periodicity in 0, but we know that 
the flux free energy is 27r periodic in 0, the external flux we add, as already visible 
in Eq. (fT^ . Therefore we have to consider the influence of all the other flux sectors, 
according to Eq. (fTTj) . and write 



- F(0) 




(46) 



F(0)-F(O) 




(47) 




(48) 



k 
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which correctly displays the periodic structure of the flux free energy. In Fig. I15lwe 
show a cartoon comparison of equations (jTrf) and (jlH)) . 




Figure 15: Qualitative behavior of the ansatz Eq. (jiSj) and Eq. (|T7j): the difference 
disappears in the limit (3 —>■ oo. 

Therefore we define Pr{P) implicitly through the equation 

^^^1^ U.o= horn ^ hi(3) (49) 

where F{(f),l3) is given in Eq. ()48p . In order to give an idea of how the map (3r{(3) 
works, observe that in the limit 13 ^ oo the following result holds 

f3n{f3) = (3[l-^ + O(-l)] (50) 

because as /3 ^ oo the renormalized coupling is equal to the helicity modulus up to 
exponentially small corrections, hence we can use the weak coupling expansion of 
SecO 

In section |S] we already showed that the ansatz Eq. (jlSj) provides an accurate 
description of the numerical data. Here we make a more refined check: this ansatz 
has to provide the right mapping between the helicity modulus at flux = and 
(f) = 71. Since it depends on only one parameter, once the curvature at the origin is 
fixed, the curvature at = vr has to be fixed too. In Table Q we check the validity 
of this assumption. As we can see, the ansatz works very well over a wide range 
of (3 values. We also observe that from the theoretical point of view, defining the 
helicity modulus as the curvature of the flux free energy at flux vr, instead of at flux 



24 
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7 / n\ 






derived from col. 2 


measured 


1.02 


0.524(1) 


-2.18(1) 


-2.16(1) 


1.10 


0.713(1) 


-4.30(1) 


-4.30(1) 


1.20 


0.852(1) 


-6.31(1) 


-6.32(1) 


1.30 


0.972(1) 


-8.35(2) 


-8.32(1) 



Table 1: Check of the vahdity of the ansatz Eq. (|48p . The first column is the bare 
(inverse) coupling; the second column is the measured value of the helicity modulus 
at flux 0; the third and the fourth columns correspond, respectively, to the analytic 
prediction of the helicity modulus at flux vr, derived from that at flux according 
to Eq. (j49|) ■ and to the directly measured quantity (the measurements are made on 
a 10^ lattice). 



0, could have been a perfectly consistent choice: 



KW) (51) 



Therefore, Eq. (j3T|) provides an equivalent way to compute the renormalized cou- 
pling. Moreover, the data obtained at flux vr are statistically more accurate. 

We measured the two helicity moduli (at flux and n) at (3 = (3^ , in the 
Coulomb phase, on a 32"^ lattice, obtaining /io,+ = 0.4325(15) and = —1.413(16) 
respectively. Then, we extracted the renormalized coupling from Eqs. (|49|) and (|5ip. 
This gives 

el = ^ = 2.30(1) for = (52) 
el = -^ = 2.31(1) for0 = 7r (53) 

PR 

which are in agreement (on smaller lattices the renormalized coupling at flux zero 
is systematically smaller than the one at flux tt; one could use this difference as 
a criterion to evaluate how close to the thermodynamic limit we are.). We take 
as our final value the average, with a conservative error sufficient to cover both 
determinations: 

el = 2.305(15) (54) 

This value can be compared with those in the literature (up to a Air factor) 
ej^ = 2.08(14) 1221 and = 2.39(12) 123- Our value is more accurate. The main 
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reason is that the other authors compute e|. from a fit of the static potential between 
two charges, obtained from Wilson loops measurements; the limitations in the size 
of the loops makes the estimate more difficult. In our case, instead (see Eq. (j2ni)), 
we do not have this limitation and use the lattice in its whole extension. 

We can now go back and consider the analysis performed in the caption of Fig. |21 
The parabolic curve, computed at flux = ±2n, gives the classical suppression of 
the corresponding flux sectors. But a renormalized (3 appears in front of it: the value 
PeS. = 0.85 was computed through the requirement that the parabola had to pass 
through the local minima. Instead, from a measurement of the helicity modulus at 
(3 = 1.2 one flnds jS^i = 0.852(2). This perfect consistency is yet another check of 
the vahdity of Eas. (l47l48ll . 

7.3 Universality of the helicity modulus 

In j2j Cardy made a conjecture about the universality of the helicity modulus (that is, 
of the renormalized coupling) at the transition point. In the typical renormalization 
group language, 'universality' means that we are going to flnd the same value in all 
Abelian theories that share the same 'long range' properties, when approaching the 
boundary with the conflned phase from the Coulomb side. 

This conjecture is motivated by the formal analogies existing between the 4d 
compact f/(l) theory and the 2d XY model. Let us start by recalling the hamiltonian 
of the lattice 2d XY model 



where p is the coupling parameter (which has dimension of energy), 6{ri) is the spin 
angle at site rj, and the sum is over flrst neighbors. This leads, in the continuum, 
to 



and in this context p is called the helicity modulus (in analogy with what we did 
in the previous section, one can argue that there is a difference between the bare p, 
which is just a coupling constant, and the renormalized pn, which takes into account 
the fluctuations of the system). 

The partition functions of both theories admit a decomposition of the form 




(55) 




(56) 



Z 



(57) 
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where Zp is the partition function of a gas of massless excitations (spin waves for 
the XY model and photons - hence the subscript P - for t/(l)) and of interacting 
- hence the subscript / - topological excitations (vortices for the XY model and 
monopole current loops for t/(l)) interacting via a Coulombic potential. 

Cardy then develops an analogy between the critical exponent 77 of the Id XY 
model, characterizing the falloff of the spin-spin correlator (in the massless, large-/? 
phase) 

and the renormalized coupling for the Abelian theory: both quantities arise from the 
additive renormalization of the bare coupling by the susceptibility of the topological 
defects 

27rr/ = ^ - 7r2^(r/a)2(m(0)m(r))/ (59) 



/5 
1 

~ ^5 24 



r 



^ e| = l-^^(r/a)2(m^(0)m^(r)), (60) 



r,fj. 



where m represents the defect field, /3 is the bare coupling and the index / in (■)/ 
represents the interacting part Zj of the partition function in Eq. (jHTjl . The crucial 
step is to show the universality of rj in the XY model. Computing explicitly the left 
side of Eq. (fSHj) [7] , one finds that 

where A;^ is the Boltzmann constant, which provides a relation between 77 and the 
helicity modulus. Taking the limit T we find ?7c = |: this result is independent 

of the specific RG trajectory considered, and therefore is universal, in the sense 
that it is unchanged when the short range features of the theory are modified. The 
combination of this result with the similarity between Eqs. (jS^ and ()60p is the basis 
for Cardy's conjecture that is a universal quantity for the Ad Abelian theory. 

The check of this hypothesis can be performed on different levels. To get a first 
qualitative idea of the analogy between the two theories we plot the quantity of 
the Ad theory fPig. I16|). which can be compared directly with the helicity modulus 
in the XY model (our definition of the helicity modulus for U{1) includes an extra 
factor P). In the XY model, the behavior near T = is given by (1 — ^T), and 
approaching by a power-law singularity {hc{l + |(1 — + 0{1 — ^)), u = 0.5) 

IZI 
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Figure 16: Helicity modulus in the 4d U{1) theory (6^ lattice), divided by P for 
direct comparison with the helicity modulus in the 2d XY model, as a function of 
"temperature" 1/(3. The linear and the power law expansions are represented around 
the regions ~ and, respectively, approaching the transition. The vertical line 
represents the transition point. 



Considering the U{1) case, we obtain the same behavior = 1 — ^ + 0{^) 
from a weak coupling expansion ^. Secondly, we measure the exponent u from the 
ansatz 

M).^(i,„(i_|n (62) 

close to the transition. The resulting exponent grows systematically with the lattice 
size, as shown in Tab. |21 From a full jackknife analysis of the data on our largest 
lattice {L = 12), we obtain u = 0.34(3). This value compares well with the value 
u = 0.43 ± 0.10 in [221, computed using the same ansatz for the renormalized cou- 
pling data (extracted from Wilson loops) in the vicinity of Pc- It is not excluded 
that the value 0.5 is attained, as in the XY model, in the thermodynamic limit. The 
overall agreement between the helicity modulus in the two theories is remarkable, 
and poses the issue of the universality of the jump of the order parameter at the 
transition on a more intuitive basis. 

^Consider the helicity modulus for the Wilson action h{P) = y{J2^^^^p) ~ ^ ^ p)'^) ' 

one can show |22 that the first term /3(cos0) = P — j + 0{^), while the second term turns out to 
be 0(4). A generalization to the extended action is straightforward. 
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L 


V 


8 


0.28(2) 


10 


0.30(1) 


12 


0.34(3) 



Table 2: Fitted exponent v from Eq. ()62j). for increasing lattice size L. 
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Figure 17: Lines of constant physics (i.e. constant renormalized Coulomb coupling 
= 1/e^) in the extended (/?,7) bare coupling plane. The numbers indicate the 
value of the helicity modulus (~ fi^) along the various iso-lines. Solid blue hues 
indicate the phase boundaries. In the confined phase surrounding the origin, the 
renormalized is zero. 

Going to a more quantitative footing, a first possible test is to compare the 
values of the renormalized coupling one obtains using different lattice discretizations 
of the continuum Abelian theory. In [22] a comparison is performed between the 
Wilson and the Villain action, which shows agreement of e|.(/3^) within errors, thus 
supporting Cardy's conjecture. 

Here we propose a more demanding test: if universality holds, the value of the 
helicity modulus should be the same along the phase boundary in the extended phase 
plane. This statement can be motivated by observing that each phase is actually 
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characterized by some specific infra-red (IR) behavior, and this is all we need to 
apply the universality hypothesis. Therefore, in order to verify Cardy's conjecture, 
we studied the lines of constant physics in the Coulomb phase of our model, i.e., 
the lines of equal renormalized coupling. They are shown in Fig. [T7I 

From this figure we can deduce that Cardy's conjecture is not satisfied. On one 
hand, in the lower right corner of the phase diagram the lines of constant physics 
tend to be parallel to the phase boundary, so numerically we should find that the 
conjecture is roughly verified (in fact a numerical test performed at /3 = 2.5 fFig fTT|) 
gives a value of the helicity modulus h = 0.37(4) in agreement with the Wilson axis 
value Eq. (jSHl))- On the other hand, if we look at values of 7 > 0, it is clear that 
the helicity modulus grows, thus disproving the conjecture. In order to exhibit this 
effect, in Fig. HSl we show the helicity modulus as a function of /5, for 7 = 0.9: it 
takes value ~ 2 near the transition, totally incompatible with the Wilson axis value 
of Section rm 
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Figure 18: Helicity modulus versus /5, at 7 = 0.9 (6'^,8'*, 10^ lattice). 

As a further comment on Fig. El it is interesting to observe that the iso-lines, 
away from the boundaries, are roughly straight with slope (in Fig. El this is 
realized in the upper right corner): this slope is easily understood considering the 
helicity modulus h{P, 7) at weak coupling, where given a value h{(3, 7) = /i we find 

p + 4:-f = h (63) 
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8 Making monopoles with fluxes 



Since we have been applying external electromagnetic fluxes to our system, we con- 
sider the natural possibility of enforcing the presence of a monopole in the lattice 
through its flux. Such a possibility is provided by adding an external flux (just 
like we did in Section ^ but using different boundary conditions, namely spatial 
C-periodic b.c.^^: when we cross the boundary of the lattice moving in a spatial 
direction we enter another copy of the lattice whose links are all conjugated 

U^ix + Liii) = U;{x), t = 1,2,3, /i = l,...,4 (64) 

The conjugation produces the inversion of the orientation of the plaquettes, and also 
the inversion of the direction of magnetic fluxes. If we try to vary the flux in such a 
lattice, and we insist on requiring translation invariance, we are faced with a severe 
restriction: because of the inversion of the flux direction, on the boundaries we need 
to impose the consistency condition 

(f) = —(j) mod 2ti (65) 

so only a flux = or tt mod 2ti can be added. In the latter case the total flux 
going out of the lattice is 2tt, and therefore the presence of one magnetic monopole 
is enforced in the lattice. What we called 'flux free energy' before can now be 
rephrased as monopole free energy, and we set 

F(0 = TT, /5) - F(0 = 0, /?) = L4F™p. (/?) (66) 

i.e., -Fmonop. is the free energy per time-slice, or, in more suggestive terms, the free 
energy per unit of monopole current. 

In Fig. El we present the numerical determination of this quantity along the 
Wilson axis. First of all, observe that it is an order parameter, as expected: the 
monopole free energy is zero in the conflned phase and different from zero in the 
Coulomb phase. The flnite size effects are clearly visible, as for the helicity modulus. 
In the Coulomb phase, as we derived in the previous section, we have flnite size effects 
proportional to 1/L (at leading order). The main qualitative difference w.r.t. the 
helicity modulus is that these lines intersect each other (at (3 ~ 1.0125). We will 
comment about this during the data analysis. 

We now make the same kind of analysis we performed for the helicity modulus: 
flrst we assume that the transition is flrst-order, therefore we flt the data of the 
monopole free energy in the range (3 < (3c, according to Eq. (j31|) . The results are: 
m+ = 0.08(2), A/V6 = 0.028(5), logX = 0.50(8), /^^ = 1.0115(3), ai ~ 14.7,^2 ^ 
(5.7)^. In Fig. 1201 on the left, we show the result of the flt: the agreement is excellent, 
and we can also observe the approximately straight behavior (in log-scale) of the 
order parameter away from j3c, signaling the flrst-order nature of the transition. 
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p 

Figure 19: Monopole free energy as a function of j3 for different lattice sizes (L = 
6,8,10). 

Going to the second-order analysis, we observe that the data show good collapse 
for 13 < [3c only (see Fig.|2nion the right) when (3 is rescaled according to Eq. (jl^ for 
a 'test' value of ~ 0.45. We also observe that the crossing of the curves of an order 
parameter at the transition is most unusual for a second-order phase transition: 
the usual behavior is that of a family of non-intersecting curves approaching the 
thermodynamic limit. 

Also in this case, therefore, we observe an indication in favor of the first-order 
nature of the transition. 

As a numerical algorithm, the chain rule (Eq. (|23|)) was used also in this case. It 
is interesting to note that all the individual ratios violate the boundary consistency 
condition Eq. we mentioned at the beginning; only the product of all of them 
(corresponding to a flux of vr) respects it. To obtain the free energy over a range of 
(3 values, we simply apply the Ferrenberg-Swendsen reweighting algorithm to each 
ratio of the chain and then perform the multiplication. 

We conclude this section with some notes about other strategies used in the 
literature to enforce the presence of a monopole in the lattice. 

The Pisa group [27] builds a magnetically charged operator which satisfies the 
following equation 

/i(y,t) |^(x,t)) = \A{x,t) + -h{x - y)) (67) 

e 
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Figure 20: Comparison of the first (left) and second (right) order scenarios for the 
monopole free energy data. 



The actual disorder parameter is built considering a discretized version of Eq. (jUHj) . 
Then the two-point function 



describes the propagation of a monopole in time from to t. Finally, considering 
the usual ansatz 



one focuses on the measurement of (yu) as a function of which should provide the 
signature for monopole condensation. 

The boundary conditions used are pure periodic in the spatial directions and 
C-periodic in time direction. Given a time-slice, the total fiux going out from the 
boundaries must be zero. This means that when a monopole is added by hand at 
position an anti-monopole is also automatically created somewhere. In our case 
it is possible to enforce the presence of only one monopole in the lattice because of 
the C— periodic spatial boundary conditions. A second remark attains the property 
of translation invariance: in our case no input about the location of the monopole 
in the system is given; this choice is taken dynamically by the system. 

In order now to compare quantitatively the two constructions, we show in Fig. |^ 
the derivative w.r.t. (3 of the monopole free energy, which is the quantity that 




D(t) = (/i(x,t)/x(f,0)) 



(69) 



D{t) = ce-*^* + 



(70) 
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directly compares to the Pisa group data. It is interesting to observe that, for lattices 
of comparable size, their peak is roughly twice as high than ours, which is compatible 
with the presence of a pair of monopoles instead of only one. Furthermore, the width 
of the peak is also much larger, probably as an effect of the discretization of the 
monopole field. 
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Figure 21: Comparison of our data (the first derivative w.r.t. (3 of the monopole 
free energy shown in Fig. IT^ with Di Giacomo et al. [77j . 

We can also compare the value of the monopole free energy in the classical limit. 
In j2Hl a computation of the classical monopole free energy in the weak coupling 
approximation is performed. We numerically computed the same quantity by just 
cooling the system (limit /5 — > cxd, obtained by minimizing iteratively the action 
while visiting the links) and measuring the classical monopole free energy per time- 
slice Scl/L^: this is the relevant quantity we compare. Our results are shown in 
Fig. E21 foi' the Wilson action. We obtain 

1 2(1 

SJU = 4.695 - -J- (71) 

Ref.j2Z|, Eq. (58), provides an expression corresponding to our {—Sc\) in the weak 
coupling approximation. Since the lattice there is of size L'^ ■ 2L, we consider 

S^JL, = -<-^-°^f + ^-^^^' = 2.525 - ^ (72) 

The rather large difference in the constant terms between Eq. fl72|) and Eq. (fTTj) 
should come as no surprise. We assign it to the discretization of the monopole field 
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Figure 22: Lattice measurement of the classical monopole free energy. 

used in the Pisa group computation: since the monopole energy is UV-divergent, 
small changes in the monopole field at short distance will have a large impact. Inter- 
estingly, the 1/L terms, which are insensitive to the monopole field discretization, 
are very nearly in the ratio 1:2, again in agreement with our interpretation of the 
two constructions. 

Jersak et. al. perform an exact duality transformation of the Villain action, 
imposing spatial C— periodic boundary conditions in the spatial direction and pure 
periodic in the time direction. Thus their system is physically identical to ours. 
They compute the two-point monopole correlation functions in the time direction 
D(t) and use the ansatz 



to extract the monopole mass M and monopole condensate (/x). They measure the 
monopole condensate as a function of (3 in the confined region and the monopole 
mass (again vs. (3) in the Coulomb phase. Assuming a second order scenario for the 
phase transition, a scaling behavior is observed for both quantities as they approach 
the (conjectured) critical point. First of all, it is interesting to observe that in Fig. QUI 
(on the right) we find a 'test' critical exponent v = 0.45 which compares very well 
with the value they find u = 0.49(4) for the critical behavior of the monopole mass 
in the Coulomb phase (Eq. (1.2) in |2H])- Moreover, from Fig. 1201 (on the left) we 
deduce a value of the monopole free energy at the transition of about 0.1, which 
is smaller than the smallest value ~ 0.3 computed in [25, just above the 'critical' 




(73) 
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point. So Ref. j2Sl was not in a position to see the jump of the monopole mass. 

9 Conclusions and perspectives 

In this paper, we have followed 't Hooft's viewpoint that twist (or flux) free energies 
can serve as order parameters for confinement. Twisted boundary conditions were 
introduced for non-Abelian SU{N) gauge theories. There, the amount of twist is 
quantized according to the center group Zjy. Here, we have considered the Abelian 
f/(l) theory, and the twist can be continuously varied. Twisted b.c. are equivalent 
to periodic b.c. plus an external electromagnetic flux; we measure the free energy of 
such a flux and show that it is an order parameter, which vanishes in the confined 
phase just like in the Yang- Mills case [23 HZl EDI- 

Studying the U{1) theory from the perspective of fluxes allows us to make prac- 
tical and conceptual progress. 

On the practical side, we have introduced the helicity modulus for a gauge theory, 
which represents the response of the system to an infinitesimal external flux. This 
observable is an order parameter just like the flux free energy. It is simple to 
measure, and gives a clear signal for the phase transition even on a small lattice, 
because it makes use of the full extent of the lattice. Moreover, it is simply related 
to the renormalized coupling. Using this observable, we show that the transition 
is first-order, with parameters {jSc, latent heat, phase asymmetry) in full agreement 
with the recent results of Note that a jump in a non-local order parameter 

is not, by itself, the signature of a first-order transition. For example, the helicity 
modulus of the XY model has a (universal) jump at T^, although the transition is 
continuous. The first-order signature comes from the exponential (not power-law) 
behavior of the finite-size effects, as seen in Fig. left. 

From the helicity modulus we extract the renormalized charge = 2.305(15) 
at the transition point with high precision. We map lines of constant renormalized 
charge in the (/3, 7) coupling plane, and disprove the conjecture by Cardy .21 of the 
universality of along the transition line. 

On the conceptual side, we have shown that a consistent description of the U{1) 
phase transition can be given without refering to magnetic monopoles, using the 
notion of flux only. In fact, we can build a monopole from flux, and force a monopole 
into the system by adjusting the boundary conditions to produce a 27i outgoing flux, 
preserving gauge and translation invariance. We compare this elegant approach with 
earlier constructions. 

The vanishing of the monopole free energy in the confined phase is thus accom- 
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panied by the vanishing of the flux free energy. Flux sectors, characteristic of the 
Coulomb phase, disappear. The condensation of monopoles is accompanied by the 
condensation of vortices or Dirac sheets. 

From our study, it is clear that monopoles are not the only choice of topological 
defect to describe the phase transition. Fluxes and monopoles are to each other like 
"chicken and egg", and assigning preeminence to one or the other is a subjective 
matter in the U{1) case. In the Yang-Mills case, fluxes are more commonly called 
center vortices, and monopoles become Abelian monopoles after Abelian gauge- 
fixing and projection. Again they are similarly related to each other. This time, 
the center vortices may be given some preference, because their free energy is a 
gauge- invariant, UV- regular order parameter [T71 1^. 

Several extensions of our work come to mind. The dimension of the system 
could be reduced to 3. The helicity modulus is still an order parameter, but at 
zero temperature it is always zero, since the theory is always confining [3T]. At 
finite temperature however, a deconfinement transition occurs |32J. This transition 
is difficult to study numerically, because according to the Svetitsky-Yaffe conjecture 
it belongs to the Kosterlitz-Thouless universality class and is of infinite order. 
The U{1) helicity modulus, which in this case directly maps onto the XY helicity 
modulus, would be a very appropriate order parameter. Another direction would 
be to measure the helicity modulus in a Yang-Mills theory after Abelian projection. 
It may give information about the deconfinement transition as the temperature is 
varied 

Finally, the picture in terms of fiuxes is the natural environment to introduce 
a U{1) topological charge. Let ^^i, be the fiux through the orientations /z, z/, and 
define $oi = ki, and e^k^jk = rrii. The topological charge Q can be expressed as 

Q = k-m (74) 

Even though instantons do not exist in the U{1) theory on i?^, the topological charge 
above, defined as the integrated second Chern character, can be non-zero on a torus 
in the presence of fiuxes.^ This would provide a natural laboratory to study the 
topological charge susceptibility or 6*- vacua j35|, at least in the Coulomb phase. 
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APPENDIX: Summary of the simulations performed 
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Table 3: Summary of the numerical simulations performed for this paper. The four 
columns contain respectively the lattice size L, the Wilson and extended couplings 
{j3, 7) (when two values of /3 or 7 are indicated, e.g. (0.99 -7- 1.02, 0), a grid of points 
is meant in between), the number of measurements, the updating algorithm ('h.b.' 
indicates heat-bath, 'o.r.' over- relaxation; for the 'snake update' [T7j, see Sec|Hl), the 
main observable. 
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